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ON  THE  METHOD  OF  GREEN'S  FUNCTION  IN 
THE  THERMOELASTIC  THEORY  OF  SHALLOW  SHELLS 


By 

R.  P.  Nordgren 


Abstract 


The  metlllod  of  Green's  function  is  applied  to  the  quasi-static 
thermoelastic  theory  of  shallow  shells  with  heat  conduction  equations 
included.  Solution  formulas  are  derived  for  middle-surface  displace¬ 
ments  and  stress  and  temperature  resultants  in  terms  of  initial  and 
edge  temperatures,  internal  heat  sources,  ambient  temperatures  at  the 
upper  and  lower  surfaces,  and  surface  tractions.  Equations  are  given 
for  the  Green's  functions  appearing  in  the  solution  formulas.  Exten¬ 
sion  to  more  general  shell  theory  is  discussed.  By  way  of  example, 
the  method  is  applied  to  thermoelastic  problems  for  two  classes  of 
shallow  shells.  Also,  the  effect  of  transverse  shear  deformatioTi  is 
examined  with  reference  to  a  shallow  spherical  shell. 

1.  INTRODUCTION 

This  paper  applies  the  method  of  Green's  function  to  quasi-static 
thermoelastic  problems  in  the  theory  of  shallov;  shells.  Thermoelastic 
equations  for  shallow  shells  follow  from  the  original  vork  of  Marguerre 

upon  addition  of  the  effect  of  thermal  expansion.  The  two  temper¬ 
ature  resultants  appearing  in  these  equations  are  governed  by  two-dimen¬ 
sional  heat  conduction  equations  derived  by  Bolotin  (i960)  for  thin  shells. 
For  flat  plates  Bolotin's  equations  reduce  to  those  of  Marguerre  (1935) 


*  The  nonlinear  terms  contained  in  Marguerre 's  equations  will  not  be 
included  here. 


and  for  shallow  shells  they  are  Independent  of  the  form  of  the  shell  middle 
surface. 

Employing  the  Indicated  equations,  solution  formulas  are  derived  which 
express  the  middle- surface  displacements  and  the  stress  and  temperature  resul¬ 
tants  in  terms  of  the  initial  temperature,  edge  temperature,  heat  sources 
within  the  shell,  and  ambient  temperatures  at  the  upper  and  lower  surfaces. 

The  two  sets  of  Green's  functions  in  these  integral  solution  formulas  are 
shown  to  satisfy  the  combined  thermoelastic  equations  with  the  two  resultant 
heat-supply  terms  replaced  by  &-functions.  A  related  solution  formula  for 
shallow  shells  under  surface  and  edge  traction  is  also  obtained,  together 
with  equations  for  the  Green's  functions  appearing  in  it. 

The  solution  formulas  obtained  extend  immediately  to  a  simplified  the(pry 
of  shells  (sometimes  referred  to  as  the  "technical"  theory  of  shells)  pre¬ 
sented,  e.g.,  by  Green  and  Zerna  (195*+)-  Extension  to  more  general  shell 
theory  is  possible  and  Is  discussed.  The  method  of  derivation  of  the  thermo¬ 
elastic  solution  formulas  is  applicable  to  other  special  theories  of  elasti¬ 
city  as  well  as  to  the  three-dimensional  theory.  In  this  connection  it 
should  be  recalled  that  Parkus  (1959»  P*  Ij)  gives  a  solutioii  formula  for 
stress  in  quasi-static  thermoelastic  problems  on  the  basis  of  analogy  and  an 
intuitive  argument.  Also,  integral  formulas  for  several  thermoelastic  prob¬ 
lems  with  unspecified  temperature  distribution  are  derived  by  Goodier  (I958) 
and  Goodier  and  Nevill  (I96I),  Sliigular  solutions  for  shallow  spherical  and 
cylindrical  shells  with  concentrated  temperature  resultants  are  given  by 
FlUgge  and  Conrad  ( I956,  I958).** 

The  solution  formulas  derived  here  are  applied  to  thermoelastic  prob¬ 
lems  £01  two  classes  of  shallow  shells,  namely,  unlimited  shallow  shells 
with  quadratic  middle  surfaces  and  rectangular  shallow  shells  with  simply 
supported  edges  parallel  to  the  principal  axis  of  quadratic  middle  surfaces. 

With  the  aid  of  a  method  of  spectral  representation  established  by  Friedman 
(1956),  the  thermoelastic  Green's  functions  are  obtained  for  the  former  class 
in  the  form  of  a  Fourier  series  of  Hankel  transform  integrals  and  for  the 
latter  class  in  the  form  of  convergent  double  Fourier  series.  Both  representations 


Previously,  one  of  these  singular  solutions  for  a  shallow  cylindrical 
shell  was  proposed  by  Rabotnov  ( 1946)  as  the  solution  for  a  concen¬ 
trated  normal  force.  However,  an  algebraic  mistake  occurred  and,  as 
shown  by  Flflgge  and  Conrad  ( I958) ,  the  concentrated  force  solution 
cannot  be  represented  by  a  singularity  of  this  type. 
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appear  suited  to  numerical  evaluation.  For  three  members  of  the  class  of 
unlimited  shallow  shells,  the  Green's  function  is  evaluated  at  the  source  point 
as  a  function  of  time  in  terms  of  tabulated  functions  and  comparison  is  made 
with  the  dynamic  thermoelastic  Green's  functions  for  bending  of  a  flat  plate. 

An  example  of  edge  heating  is  studied  for  the  class  of  rectangular  shallow  shells. 

In  the  last  section,  the  effect  of  transverse  shear  deformation  is 
examined  with  reference  to  an  unlimited  shallow  spherical  shell  (parabaloid 
of  revolution)  under  specified  temperature  field  and  normal  surface  traction, 

A  quasi-static  solution  is  obtained  by  the  method  of  Green's  function  with 
use  of  the  extended  definition  of  the  Laplacian  operator  following  Friedman 
(1956)*  The  character  of  the  fundamental  singularities  is  compared  with  that 
of  previous  solutions  by  Reissner  (I9I46)  and  by  FlHgge  and  Conrad  (195^)  according 
to  classical  theory,  which  neglects  the  effect  of  transverse  shear  deformation. 

2,  QUASI -STATIC  THERMOELASTIC  THEORY  OF  SHELLS 


2.1  A  Solution  Formula  for  Heat  Conduction.  The  two-dimensional  heat  con¬ 
duction  equations  for  thin  shells  derived  recently  by  Bolotin  ( I96O) 
from  a  general  variational  principle  for  three-dimensional  Isotropic 
(uncoupled)  heat  conduction  with  linear  heat  transfer  at  the  boundary 
may  be  written  as 


0  .  1  ^ 
N  'n  N  K  bt 

-  n..e .  _  i  ^ 


(2.1) 
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In  (2.1)  the  Laplaclan  is  V  V„  ,  where  the  covariant  derivative 

Q*  P  Q, 

operator  V  refers  to  a  curvilinear  coordinate  system  x  (a  =  1,  2)  on  the  middle 
surface  of  the  shell  with  metric  a  and  conjugate  metric  ,  the  summation 

Up 

convention  being  understood.  The  thickness  of  the  shell  is  denoted  by  h  and  x® 
is  the  coordinate  normal  to  the  middle  surface.  The  temperature  resultants  6^^ 
and  6^  are  expressed  in  (2.2)  as  Integrals  of  the  temperature  T  (above  a  fixed 
reference  temperature),  and  the  derivation  by  Bolotin  (i960)  is  based  on  the 
assumption 

T  »  0j,(x“)  +  x2e^(x°').  (2.5) 

The  heat  inputs  Q  and  Q  arise  either  from  internal  heat  sources  of  rate  Q 
H  N 

per  unit  time  per  unit  volume  or  from  heat  flux  across  the  upper  and  lower  sur¬ 
faces  of  the  shell  where  the  ambient  temperatures  are  and  respectively. 

The  surface  conductance  H,  thermal  conductivity  K,  specific  heat  c,  and  mass 
density  p  are  assumed  independent  of  temperature  as  well  as  coordinates.  On  the 
boundary  curve  C  of  the  shell  middle  surface  the  temperature  resultants  may  be 
specified,  i.e., 

which  is  a  special  case  of  the  linear  heat  transfer  edge  conditions  given  by 
Bolotin  (i960).  In  addition,  and  must  satisfy  the  initial  conditions 


lim  0jj(x  ,t)  =  0“(x“), 
t  -♦  o 


lim  0„(x“,t)  =  e^(x“). 


(2.5) 


where  T°  «  0°  +  is  the  Initial  temperature  of  the  shell. 

Since  the  heat  conduction  equations  (2.1),  (2.4),  and  (2,5)  are  of  the 

same  form  for  6.,  and  6^,  we  need  only  discuss  solutions  for  ;  solutions  for 
N  H  N 

then  follow  merely  by  replacing  the  subscript  N  by  M,  We  shall  use  the  linear 
vector  space  of  all  functions  such  that 


0  dS  <  CO  , 

N 


(2.6) 


and  the  scalar  product  of  two  functions  and  0^  in  this  space  is  defined  as 


(2.7) 


where  dS  denotes  the  element  of  surface  area  on  the  shell  middle  surface. 


Equation  (2»l)  contains  the  operator 


T  .  i  ^ 

S  "  ^  ~  \  “  K  bt 


(2.8)' 


whose  domain  will  be  the  set  of  all  functions  9  in  the  vector  space  such  that 

N 

Sx  and  /St  are  piecewise  continuous  and  of  Integrable  square  and 

N  Q!  p  N 

such  that  6  satisfies  the  zero  Initial  condition  lim  9  =  0  and  the  homo- 

N  t  -»  o*^ 

geneous  boundary  condition  0^^  =>  0  on  C.  With  the 

theorem  for  surfaces,  we  may  establish  the  identity 

lyi  - 


aid  of  Green's 


iH 

■U 


(2.9) 


1  im  9  9  “1  im  9  9 

Aim  lim 


dS, 


cc 

where  n  is  the  normal  to  C  and 


(2.10) 


is  the  operator  adjoint  to  whose  domain  is  the  same  as  that  of  L^,  except 
that  instead  of  a  zero  initial  condition,  lim  0, 


t!  ->  00 


N 


N’ 

0.  Thus,  the  right-hand 


side  of  (2.10)  vanishes  when  0^^  and  0^^  are  in  the  domains  of  and 
respectively.  The  definition  of  Lj^  may  be  extended  to  functions  0|J  not  in  the 


domain  of  L.,  by  letting 


(2.11) 


where  1.0"  may  be  a  symbolic  function. 

N  N 

We  shall  show  that  the  Green's  function  for  L„  is  the  actual  or  symbolic 

M 

function  Gj^(x,  t;x',t')  which  satisfies 


LjjGjj(x,t;x',t')  -  -  6(x  -  x')  5(  t  -  t')  , 


Gj,(x.t;x',t')  =  0 
Gj,(x,t;x',t')  =  0 


if  t  <  t'  , 
X  on  C  , 


(2.12a) 

(2.12b) 

(2.12c) 


where  6(x,-  x')  is  the- two-dimensional  B-function  and  8(t  -  t')  the  one- 
dimensional  B-function. 
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Replacement  of  t  by  -t  in  (2.12)  results  in 


LjjGj^(x,  tjx",  t")  =  -  6(x-x")8(  t-t") 
G*(x,t;x",t")  =0  if  t  g  t", 

G^(x, t;x", t")  =0,  X  on  C, 


(2.13a) 

(2.15b) 

(2.13c) 


where 


Gjj(x,t;x",t‘')  =  Gjj(x,-t;x",t") 


(2.U) 


can  be  shown  to  be  the  Green's  function  for  the  adjoint  operator  L  .  Identi- 
fication  of  with  Gjj(x,  t;x' ,  t ')  and  0^^  with  Gj^(x,  t;x",  t")  in  (2.9)  leads  to 
the  relation 


Gjj(x",t";x',t')  =  G^(x'.t';x",t"). 


(2.15) 


while  identification  of  0'  with  G  (x, t;x',t')  and  use  of  (2.I5)  results  in 

N  N 

the  following  solution  formula  for  0^^  in  terms  of  its  Initial  values,  boundary 
values,  and  Oj^i 

ej^(x,t)  dt'  G^(x,  t;x',t')  Qj^(x’,t’)dS' 

-  f  dt'  r  n“(x')V  gJx,  t;x',t')0  (x',t')ds'  (2.16) 

J  Q  .)  Q,  aw  n 

+  I'  ^  G^(x,t;x',0)  0j^(x',O)dS'. 

A  formula  for  0^^  follows  by  replacing  subscript  N  by  M  in  (2.I6). 


2.2  A  Solution  Formula  for  Shallow  Shells.  We  recall  that  the  stress  differ¬ 
ential  equations  of  equilibrium  for  shallow  shells  may  be  written  as 


VpiPP  +  =  0. 

Vb  =0  , 


Q  =  0, 


(2.17a) 

(2.17b) 

(2.17c) 


where  q“.  p“  and  p  are  the  stress  resultants,  shear  stress  resultants, 

stress  couple  resultants,  tangential  surface  tractions,  and  normal  surface 
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traction,  respectively;  H  is  the  distance  of  the  shell  middle  surface  from  a 
reference  plane,  and  the  covariant  differentiation  now  refers  to  a  curvi¬ 
linear  coordinate  system  in  the  reference  plane.  The  strains  7  and  changes 

Cep 

in  curvature  k  of  the  middle  surface  of  the  shallow  shell  are  given  by 

0!p 


27  ^  =  V„v  +  V  -  2V  V„Zw, 
'a^  pa  a  p  a  p  ’ 


(2.18a) 


’'ap  “  "'^a^p'^’ 


(2.18b) 


where  v^  and  w  are  the  tangential  and  normal  displacements  of  the  middle  surface 
of  the  shell.*  The  constitutive  equations  for  isotropic  homogeneous  shallow 
shells,  with  the  effect  of  thermal  expansion  included  and  the  effect  of  transverse 
shear  deformation  neglected,  take  the  form 


where 


t 


Eh  ap  a 


.  Eh»  a(S  g 
12^  \i]  12(l-v  )  “  M’ 


(2.19a) 


(2.19b) 


and  Young's  modulus  E,  Poisson's  ratio v  and  the  coefficient  of  linear  thermal 
expansion  a  art  assumed  Independent  of  temperature  as  well  as  coordinates. 

Upon  substitution  of  (2.18)  and  (2.I9)  into  (2.1?)  we  obtain  the  relation 


rP  ^  hi 


i3\q  ■  I2(l-v) 


aa“Pve„ 


(2.21) 


and  the  displacement  equations  of  equilibrium 


[VA  ■  +  1’°  ■ 

-  +  V  V  v.  -  V.  V  Zw 

a  p  L  h  ^  ^  q 


-  V.  V^Zw 

T]  A  ^  ^ 


(2.22a) 


(2.22b) 


*  Marguerre  (I938)  and  Green  and  Zerna  (195^)  employ  displacements  normal  and 
parallel  to  the  reference  plane  (axial  ai;d  longitudinal  displacements)  which 
are  related  to  the  normal  and  tangential  displacements  w  and  v  used  here  and  in 
other  works  by  w  and  -  V^Zw,  respectively,  to  within  the  approximations  in¬ 
herent  in  the  theory  of  shallow  shells.  Also,  the  second  fundamental  form  of 
the  shell  middle  surface  has  the  coefficient  b  =  V^V^Z  to  within  the  accuracy 


of  shallow  shell  theory. 


ap  a  p 
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-  a'  (v^zQj^  +  +  P  =  0, 


(2.22b) 


where 


n  Eh-^ 

°  "  12(1^) 


a-  = 


XT^  • 


With  the  aid  of  (2.I7)  and(2.l8)  and  the  divergence  theorem,  it  is  easy  to 
establish  the  relation 


(2.23) 


where 


I  =  J  ■  >^*^VpW  +  q“ w)n^ds, 


(2.24) 


and  n  is  now  the  normal  to  the  curve  C  bounding  the  region  R  occupied  by  the 

(jC 

Ct 

projection  of  the  shell  middle  surface  onto  the  plane  of  the  coordinates  x  .  In 
the  classical  theory  of  shallow  shells,  since  the  effect  of  transverse  shear  de¬ 
formation  is  neglected,  the  line  Integral  in  (2.24)  should  be  written  in  an 
alternate  form  to  reflect  the  fact  that  only  four  boundary  conditions  may  be 
prescribed  on  C.  Thus,  following  Green  and  Zerna  (195^)»  write  the  components 
of  the  couple  resultant  on  C  for  directions  normal  and  tangential  to  C  as 

»» "  «r  ■ 


respectively,  and  for  future  convenience  we  do  the  same  for  the  components 
of  stress  resultant  and  tangential  displacement,  l.e.. 


"n  ■  Vs  ■ 


“t  •  ■ 

^  . 


(2,26) 


where  e  and  are  tensors  of  the  e- systems.  Denoting  the  derivatives 

pA. 

along  the  normal  to  C  by  8/9n  and  S/Ss,  respectively,  and  defining  the  perhaps 
symbolic  derivative**  through 


(2.27) 


**  See  Friedman  (I956),  pp.  140-142. 
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(2,24)  becomes 

"  'L  [Vs  -  Vt  ■  “»  S  +  ’”] 


where 


(2.28) 


(2.29) 


may  be  interpreted  as  the  total  shear  resultant  on  C.  If  Hj  has  a  jump 
discontinuity  of  magnitude  AHj,  at  s  *  s',  then  by  (2.27)  and  (2,29) 

V=Q%^+^  +  AIij.8(s-s'),  (2.50) 

where  ^  denotes  differentiation  in  the  ordinary  (not  symbolic)  sense.  The 
last  term  in  (2.50)  represents  a  concentrated  load  which  is  a  well-known 
result  in  classical  shell  (and  plate)  theory.  Appropriate  boundary  conditions 
for  shallow  shells  are  apparent  from  (2.28)  and  may  be  listed  as 


(2.31) 


either 

ea  N 

N 

or 

''n 

either 

* 

=  Nt 

or 

either 

* 

or 

Sw 

either 

V 

=  V* 

or 

w 

on  C,  where  the  starred  quantities  are  boundary  values.  In  a  later  section 
we  shall  consider  an  example  in  which 

^  ’  ‘Si  “  ’  v,j,  =  w  =  0,  on  C,  ('2.32) 

and  these  boundary  conditions  will  be  used  to  Illustrate  the  method  of  the 
present  section.  By  (2.I9),  (2.23),  and  (2.28),  the  reciprocity  theorem 
(Green's  formula)  for  shallow  shells  is 
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dS 


(2.33) 


where  primed  end  unprlmed  variables  denote  any  two  solutions  of  the  shallow 
shell  equations.  The  thermoelastic  reciprocity  theorem  (2.55)  may  also  be 
reached  from  the  isothermal  reciprocity  theorem  [0^  =  =  0  in  (2.53)]  by 

means  of  an  analogy  between  the  isothermal  equations  and  the  thermoelastic 
equations  of  shallow  shells."^  By  (2.I9)  to  (2.22),  the  analogy  may  be  written 


Pc  =  P  -  a'(v^zS„  +  12 

=  l^P+a'a“%  .  (2.34) 

+  g  a'J^\  . 


where  subscript  zero  denotes  a  variable  in  the  isothermal  equations.  As  an 
aid  in  discussing  solutions  of  (2,22)  it  is  convenient  to  denote  components 
of  the  middle  surface  displacement  vector  by  (i  =  1,  2,  5)  where 

”  '"a  •  (a  =1,2),  U3  =  w,  (2.55)'^ 


and  consider  the  linear  vector  space  of  all  such  that 


(a 


ap, 


V  v„ 
a  p 


+  w^)dS  <  00 


(2.36) 


The  scalar  product  of  two  displacement  vectors  and  in  this  vector 
space  is  defined  as 


( a^^v  v'  +  ww')dS  , 

OJ  p 


(2.37) 


"^The  analogy  is  similar  to  that  given  by  Duhamel  ( I838)  for  three-dimensional 
thermoelasticity,  and  the  Betti  reciprocity  theorem  has  been  extended  to 
thermoelasticity  with  the  aid  of  this  analogy  by  Goodler  (I958). 

4*4* 

The  only  tensor  properties  of  are  those  of  its  components  v  and  w 
vhich  are  a  vector  and  scalar,  respectively,  with  reference  to  t’^ansformations 
of  coordinates  in  the  reference  plane. 
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In  (2.22)  there  appears  the  operator  defined  by 


extended 


n 

to  functions  not  in  its  domain  by  the  definition 


dS 


dS. 


(2.44) 


We  shall  show  that  the  Green's  function  for  is  the  perhaps  symbolic 
(tensorial)  function  Gj^j(x;x‘),  which  satisfies  the  equation 


Gjk(x,x') 


■\  8(x-x')  , 


(2.45) 


and  homogeneous  boundary  conditions*  corresponding  to  (2.32)  which  with  the 
aid  of  (2.19)  and  (2.26)  may  be  written  as 

^Ni  “  ■  G3,i(x;x-)v^V^Z]  =  0  , 

“Ni  "  ^  -  0  ,  (2.46) 

»  n“  G^^(x;x')  =  0  ,  G^^(x;x')  =  0  , 

for  X  on  C.  In  (2.4l),  identification  of  Uj  with  G^j^  (x;x')  and  Uj  with 
Gj2^(x;x")  followed  by  application  of  (2.45)  and  (2.46)  yields  the  relation 

G^^(x';x")  =  Gjj^(x",x')  .  (2.47) 

Similarly,  identification  of  Uj  with  Gjj^(x;x')  in  (2.4l)  and  application  of 
(2.32),  (2.59),  (2.45),  (2.46),  and  (2.4?)  results  in 

^i(x)  P^(x')  +  Tl  ^  (®(x'))  (^,j(x;x')]  dS' 

(2.48)** 

+  [fi^(x')G^p(x;x')N*(x')  -  G^^(x;x')  H^(x')]  ds'. 


It  should  be  noted  that  under  some  boundary  conditions  in  the  original  problem, 
the  boundary  conditions  on  the  Green's  function  will  not  be  homogeneous.  For 
example,  if  and  N.j.  are  prescribed  on  C  in  the  original  problem,  corresponding 
homogeneous  boundary  conditions  on  the  Green's  function  would  not  be  possible  in 
view  of  over-all  equilibrium  requirements.  In  such  cases  the  Green's  function 
must  satisfy  inhomogeneous  boundary  conditions  which  may  conveniently  be  taken 
in  a  form  similar  to  rigid  body  displacements  as  done  for  a  flat  plate  by  Berg¬ 
man  and  Schlffer  (1953j  P*  239). 

The  primed  operators  act  on  the  x  variables  of  the  Green's  function. 
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The  solution  formula  (2.48)  expresses  the  displacements  In  terms  of  the  surface 
and  edge  tractions  and  the  temperature  resultant  field.  In  the  absence  of  sur¬ 
face  and  edge  tractions  (2.48)  may  be  written  as 


Uj^(x)  =  J  (x;x')  6j^(xl)  +  0^  (x;x')Sjj(x')]  dS', 


where  by  (2.42) 

(x;x')  =  a'  (x'>7^  G^(x;x')  -  y  (x')  Gj^3(x;x')  ]  , 

(x;x')  =  -  G^(x;x')  . 


(2.49) 


(2.50) 


By  manipulation  of  (2.45)  it  Is  easy  to  show  that  (x;x' )  satisfies  (2.22) 

with  0  =  8(x-x')  and  9  =  p^  =  P  =  0,  l.e., 

N  M 


(x;x')  =  a*  a^^y  &(x-x')j 

j  P 

Uj  (x;x')  =  ofy^z  6(x-x'), 


(2.51a) 


while  (x;x')  satisfies  (2.22)  with  &  =  6(x-x')  and  & 

1  n  N 

TpJllJ  (x;x-)  =  0, 

(x;x')  =  a'  72  5(x.x’). 


a 


P 


0,  l.e., 


(2.51b) 


Further,  by  (2.46),  satisfies  the  homogeneous  boundary  conditions 

-  u"  (»!k')  )  -  0, 

fVrvN  /  \  'vN  y  V 

"  Uq,  (x;x')  =  U3  (x;x')  =0,  X  on  C, 


(2.52) 


and  Identical  boundary  conditions  are  met  by  ^  .  The  stress  and  couple 
resultants  may  now  be  obtained  from  (2.I9)  and  (2.48)  or(2.49). 


2.3  A  Solution  Formula  for  Combined  Thermoelastic  Problems  of  Shallow  Shells. 
The  solution  formula  (2.16)  for  heat  conduction  and  the  elastostatlc  solution 
formula  (2.49)  be  combined  and  after  Interchange  of  the  order  of 
integration  written  as 
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(2.53) 


U^(x,  t)  =  j  dt'J  ^  (x,t;x|t')  Qj^{x',t') 

(x.tjx'.f )  Q„(x>,t' )]  dS'  -jT’" 

J'  n“(x')  (x,t;x',t')  0|(x',t')  +  V^^(x,  t;x' ,  t '  )e*(x' ,  t ' ) 

c 

■jT^  ji^NCx.tjx'.O)  e;(x')  (x,t;x'.0)  0»(x') 


ds ' 


dS', 


where 


(x.t;x",t")  =  f  uJJ  (x;x')  G  (x'.t;x",t") 

i  V-/  D  » 


dS' 


(2.5l^) 


(2.55) 


and  similarly  for^^  .  The  truth  of  (2.53)  easily  verified,  since  by 
(2.51a)  and  (2.5U) 

(x.t;x",t")  -  a'  a“P  (x,t;x",t"), 

(x,t;x",t")  =  a*  Z  (x,t;x",t"), 

(x.t;x",t")  =  0  , 

(x,t;x",t")  .  ^a'  (  x,  t;  x",  t”) , 

and  then  by  (2.53)  and  (2.16) 

Uj  (x,t)  =  a'a“P  Vp  (x,t), 

Uj  (x,  t)  =a'  Z  0j^(x;t)  +  ^ 


which  are  precisely  the  field  equations  (2.39)  with  P  =  0.  By  (2.54), 
and«^^  satisfy  boundary  conditions  of  the  form  (2.52).  Thus,^^  and 
play  the  role  of  Green's  function  for  combined  quasi-static  thermo¬ 
elastic  problems  in  the  theory  of  shallow  shells  and  (2.53)  determines  the 
middle  surface  displacement  in  terms  of  the  surface  or  internal  heating, 
edge  temper-nture,  and  initial  temperature.  By  (2.I6),  (2.I9),  and  (2.55)  i^he 
stress  resultants  are 
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(x,t)  =  j  J  ^  (x,t;x',t')  Qj^  (x',t') 

+  (x,t;x',t')  (x',t') 

-f  J"  n“(x')  li^^(x,t;x',t')ejj  (x',t') 

+  t;x',t')  0jj(x',t') 

+  f  ^  (x,t;x',0)  6“  (x') 


(2.57) 


+  (x,t;x',0)  ®^(x') 


dS'  , 


where 


(x,t;x',t')  =  f\^!!  (x,t;x',t') 


-a3(x,t;x',t')V^V^  Z  (x) 


-  a'a“^(x)  G^(x,t;x',t'). 


(2.58) 


and  similarly  for  •  By  (2.58),  is  just  the  resultant  obtained 
directly  from  the  Green's  functions  by  an  equation  of  the  form  (2.19a). 

Similar  expressions  may  be  obtained  for  the  couple  resultants  and  shear  stress 
resultants.  The  question  of  existence  of  solution  for  specified  edge  tem¬ 
perature,  initial  temperature,  and  heating  will  not  be  considered  here, 
although  some  restrictions  on  the  region  R  and  on  the  specified  functions 
certainly  will  be  necessary. 


2.k  Remarks  on  Extension  of  the  Results  to  General  Theory  of  Shells.  The  method 

developed  in  the  foregoing  may  also  be  applied  in  the  more  general  theory  of 

shells  which  are  not  necessarily  shallow  and  in  which  the  effect  of  transverse 

shear  deformation  may  also  be  included.  The  results  obtained  for  shallow 

shells  apply  immediately  to  the  simplified  theory  (sometimes  referred  to  as 

the  "technical”  theory)  of  shells,  presented,  e.g,,  by  Green  and  Zerna  , 

the  only  change  being  that  x  are  regarded  as  middle  surface  coordinates 

rather  than  coordinates  in  a  reference  plane.  With  only  alight  modification, 

the  formulas  of  this  section  also  extend  to  the  version  of  shell  theory 

*4* 

commonly  known  as  Love's  first  approximation,  in  which  the  constitutive 

^There  are  at  least  three  different  versions  of  Love's  approximation  in  the 
literature.  We  refer  here  to  the  one  given  in  vectorial  form  by  Reissner  (I941). 
A  full  discussion  and  additional  refer^^ces  are  given  by  Naghdi  ( I962) , 


equations  may  be  written  as 


(V,v^  -  , 


(2-59) 


where  the  covariant  differentiation  now  refers  to  a  curvilinear  coordinate 
system  on  the  shell  middle  surface  and  is  given  by  (2.20)  with  a*^^ 

now  being  the  conjugate  metric  for  the  middle  surface  coordinate  system. 
When  (2.59)  is  combined  vjith  the  general  equilibrium  equations  for  shells 


we  have 


V>.  -  V’  *  ^ '"p  [Vn " 

t  V^(b;  v^)]  -  c'  [,“e  -  f|  b“P  VpSj  +  p“  .  0  . 

-  f  [v^v^»  +  (bj;v^) 

+  a"®-"  (V^v^  -  b^^»)  -  «■  [f|  +  b“  e^J  +  P  .  0 


(2.60) 


which  are  of  the  form  (2.59)  with  the  operator  and  the  temperature 
terra  no  longer  given  by  (2.58)  and  (2.140).  Also,  it  is  easily  verified 
that  the  form  (2.41)  remains  unchanged.  Since  the  heat  conduction  equations 
(2.1)  with  solution  formula  (2.16)  are  valid  to  the  same  degree  of  approxi¬ 
mation  as  (2.59),  it  follows  that  the  result  (2.55)  holds  also  for  Love's 
first  approximation.  Hoviever,  the  Green's  function^,  (x,t;x',t')  will  now 
be  a  solution  of  (2.6o)  operated  on  by  L  with  L  ©  ■  a:6(x-x')  6(t-t'), 

CK  ri  ii  ii 

©u  =  p  =  p  =0  and  under  appropriate  boundary  conditions.  A  similar 
statement  holds  ^Iso  for  ^j^jin  (2.48).  In  connection  with 

Love's  first  approximation,  it  may  be  recalled  that  integral  formulas  for 
displacements  and  their  derivatives  for  the  isothermal  case  are  given  by 
Naghdl  (i960);  these  formulas  are  contained  in  (2.48)  and  its  derivatives. 

It  should  be  remembered,  however,  that  certain  inconsistencies  are  present 
in  the  constitutive  equations  (2.59)  in  that  and  do  not  vanish 
identically  under  rigid  body  displacement,  nor  do  they  satisfy  the  equation 
of  moment  equilibrium  about  a  normal  to  the  shell  middle  surface,  except 

*4^  u 

for  spherical  shells.  While  the  degree  of  error  is  of  the  order  —  (R  being 


The  inr.onsistencies  mentioned  are  also  present  in  the  constitutive  equations 
of  shallow  shell  theory,  although  they  are  not  considered  serious  in  view 
of  the  assumption  of  shallowness. 
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the  minimum  principal  radius  of  curvature),  and  ^  «lis  a  basic  postulate 
of  thin  shell  theory,  nevertheless  difficulties  have  been  encountered  in  some 
problems  and  at  least  on  theoretical  grounds  improved  constitutive  equations 
are  desirable.  To  obtain  improved  theories,  additional  higher  order  geo¬ 
metrical  terms  have  been  retained  in  numerous  derivations  of  constitutive 
equations  and  the  (often)  equally  important  effect  of  transverse  shear  de¬ 
formation  has  sometimes  also  been  included,  A  more  complete  discussion  and 
a  detailed  comparison  of  various  constitutive,  equations  is  given  by  Naghdi 
(1962).  When  more  exact  constitutive  equations  are  applied  to  thermoelastic 
problems  of  shells,  the  use  of  heat  conduction  equations  of  corresponding 
accuracy  seems  desirable.  Thus,  the  heat  conduction  equations  for  thin 
shells  derived  by  Bolotin  (i960)  could  be  Improved  by  retaining  geometrical 
quantities  of  higher  order  in  h/R,  in  which  case  it  may  also  be  necessary 
to  modify  the  assumption  (2.3).  Further,  the  temperature  resultants  6^^  and 
may  not  enter  the  constitutive  equations  in  such  a  simple  manner  as  in 
(2.59).*  While  additional  effort  may  be  required  to  extend  the  present 
method  to  thermoelastic  problems  in  improved  theory  of  shells,  the  main 
ideas  of  the  Green's  function  approach  will  still  apply  provided  a  recipro¬ 
city  (Green's)  formula  is  available  In  simple  form,  as  is  the  case  for  a 
system  of  improved  constitutive  equations  given  by  Naghdi  (I962),  In  fact, 
improved  equations  due  to  FlUgge  for  circular  cylindrical  shells  are  employed 
by  Goodier  and  Nevill  (I96I)  to  obtain  formulas  somewhat  analagous  to  (2»49)j 
and  they  also  derive  similar  formulas  for  various  specific  problems  in  theories 
of  thin  bars,  thin  plates,  membrane  shells,  and  three-dimensional  elasticity. 

In  these  as  well  as  other  problems,  it  should  be  possible  to  combine  the 
integral  formulas  for  displacements  in  terms  of  the  temperature  field  with 
integral  formulas  for  solution  of  appropriate  heat  conduction  equations, 
thereby  obtaining  formulas  for  displacements  in  the  combined  thermoelastic 
problem  as  in  the  foregoing  treatment  for  shallow  shells. 

3.  APPLICATION  OF  THE  QUASI-STATIC  THEORY; 

THERMOELASTIC  GREEN'S  FUNCTIONS  FOR 
TWO  CLASSES  OF  SHALLOW  SHELLS 

3,1  A  Class  of  Unlimited  Shallow  Shells.  First  we  shall  obtain  the  Green's 

functions  Gj^  and  Gj^  for  temperature  resultants  in  an  unlimited  shallow  shell, 

- - - . . 

Note,  however,  that  the  6^^  and  6^  terms  in  (2.59)  do  meet  the  equation  of 

moment  equilibrium  about  a  normal  to  the  middle  surface. 


-17- 


These  Green's  functions  are  independent  of  the  form  of  the  middle  surface, 
since  a  does  not  enter  Although  Gj^  and  could  be  obtained  from 

the  known  Green's  function  for  the  two-dimensional  heat  conduction  equation  by 
a  simple  transformation,  for  Illustrative  purposes  they  are  derived  here 
following  the  method  established  by  Friedman  (1956,  P-295)  since  this  same 
method  will  be  applied  later  to  obtain  the  thermoelastic  Green's  functions 
-dS"  ■  Thus,  we  may  write  (2.12a)  as 


Gj,  (x.t;x',t') 


-  fi 


-  1 

!  or 

“  ic  St 


8(x  -  x')  5  (t  -  t')  ,  (5.1) 


and  interpret  the  result  with  the  aid  of  appropriate  spectral  represen¬ 
tations  for  the  Laplacian  (in  the  entire  plane)  and  ^  (on  0  S  t  S  «>  with 

♦“if  ^ 

zero  initial  condition),  namely 

6(x  -  X')  r  r  .  (3.2) 

V  —00  w  — OQ 

5(1-1,)=^^  ds  ,  (5.5) 


ds  , 


(5-5) 


a  -  ioo 


respectively,  where  x^^  and  xg  arc  rectangular  Cartesian  coordinates.  By 
(3.1)  to  (3.3)  we  have 

^  r“  oa  +  ico  ^s(t-t')^i|j^(xj^-xj;)^l|a(x2-x^) 


p  00  AOO  p  a 


s  -f  /c(t1j^  +  If  +  li) 


d8d|3_d|dg 

(5.if) 


which,  upon  introducing 


_ ^ _ r 

+  '<(%  +  Si  52)  Jq 


exp  r  -  ts  -  tK{y\^  +  li  +  I2)  dt,  (5.5) 


r^  »  (xi-xi)2  +  (xg  -  x^)2  , 

after  some  elementary  manipulations,  becomes 

r  .-e 

Gj^(x.t;x',t')  =.  exp  (t-f)  ,  t  -  f  S  0, 

(5.6) 

=0,  t  -  t'  <  0, 

and  a  similar  equation  holds  for  G„  (x, t;x',t').  Next,  we  consider  the  class 
of  unlimited  shallov  shells  with  middle  surfaces  of  arbitrary  quadratic  form 
\  hich  by  an  orthogonal  transformation  of  coordinates  may  be  brought  to 


noted  by  Friedman  (I956),  (5*2)  and  (3.3)  contain  the  tvo-dimensional 
Fourier  transform  theorem  and  the  Laplace  transform  theorem,  respectively. 


2  =  i  * 


where  and  bs  are  the  principal  curvatures  of  the  middle  surface  of  the 
shallow  shell.  Tor  the  class  (j.?)*  Marguerre 's  (I938)  formulation  for  elasto- 
static  problems  of  shallow  shells  in  terms  of  w  and  an  Airy  stress  function  F| 
in  the  absence  of  surface  tractions  and  the  presence  of  temperature  resul¬ 


tants,  reads 


F  +  Eh  (bsLj^  +  bj^Lg)  w  +  Eh  ■  0  , 

D  w  -  (bgLi  +  bj^Lp)  F  +  (1-v)  D  “  0  , 


where 


T 

, 


"2  “  ^2^  ’ 


(3.8a) 


(3.8b) 


and  the  stress  resultants  are 


Nil  “  LgF  .  N22  =  LiF,  Ni2  »  - 


(3.9) 


By  (2.18),  (2.19))  (3.7),  and  (3.9)  the  tangential  components  of  middle 
surface  displacement  to  within  a  rigid  body  displacement  may  be  written  as 


r’'"  ri 

'Jo 


(L2F  -  vl*i?)  +  b^w  +  Of®, 


nJ 


(3.10) 


and  an  expression  for  V2  given  by  (3>10)  with  subscripts  1  and  2  Interchanged. 
To  obtain  the  thermoelastic  Green's  function^^^  (x,t;x',t'),  we  replace  F, 
w,  ®jj,  and  in  (3.8)  by  ,  Gjj,  and  0,  respectively,  whence,  with 

the  aid  of  (3*l)»  (3*8)  yields 


“-ft - - n - m  ' 

V®  +  k2(b2Li  +  biLa)^  ^  fj- 


(3.11) 


where 


1  2  Eh  12  /-  .i2\ 

^  “  T  “  h^  '  • 


ror  a  shallow  shell  with  middle  surface  analytic  at  the  origin,  (3.7)  re¬ 
presents  the  first  significant  term  of  the  series  expansion  of  E  in  the 
neighborhood  of  the  origin. 
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This  result  may  be  interpreted  using  (3.3)  spectral  representation 


■  ■  s/ 

^  B 


«1)  aj. 


(3.12) 


for  the  operator  L  and  a  similar  representation  for  Lg.  Thus,  following  the 
procedure  used  to  obtain  (3.6),  we  find  that 


jljN  kKd  -'<TiN(t-t')  ^  V 

«  3  “  gpj  i  Mil  -  I2)  . 


, 


(3.13) 


where 


(  Ill  ^  (eM  6§)  Si)(t-t:J 

I2J  ^.00  (|2+  ||)  +  ik(b2li2+  bill) 


(3.14) 


i|l(Xi  -  xl)  Msixs  -  x^) 


d|i  dig 


are  complex  conjugate  functions.  The  component^^i  of  the  Green's  function 
by  (3.10),  (3.4),  and  (3.I3)  is  ^ 


^  ^  (t-f)  r  “  r  “  L.  Liii(xi-x^)_  ^  iiix'  ■ 

^  'J.co  J.ce  L 

(E?  +  E§)  (li  -vlf)  +  bik^  (Ei^+  li)(bg|g  bill)  _  ^ 

(5i‘'+  li)"*  +  4^  (b^li  +  bila^)^ 


(3.15) 


bp|f  +  bi 


4(|:f  +  li)(t-t')  ^i|g(x2  -  xp  ^ 


N 

and^s  (x*b;x',t')  is  given  by  (3.I5)  with  subscripts  1  and  2  interchanged. 

.  M 

Following  the  same  procedure  as  before,  the  Green's  functionsg^  are 
found  to  be 


rM  ^  (l+vl 

1  4«' 


jca  ^-KTiM(t-t')  p+  “  r+  ”  1  ['^i|i(xi-xi)_  ^-l£ix'‘ 

'-w  SO  ^  00 


(3.16) 


(Ef  +  £i)(b2i?  +  biiDdi  -  vlf)  - 
•  (S!  +  li)"  +k2(b2|?  +  bi|2e)2 


t2j.  »5>3 


+  li)(t-t')  /|2(x2-xi)  ^ 
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(3.17) 


M 

a  similar  expression  for^  2  (x,  t}x',t'),  and 

Next,  in  order  to  study  the  integrals  and  defined  by  (3.11-)  we  intro¬ 
duce  polar  coordinates  through 


■  p  cos  ^  , 


Xj^  -  x{  =  r.  cos  6  , 


la  ■=  p  sin  4'  , 


Xg  r  sin  Q  , 


(3.18) 


and  then 


( S  fj. 


«  r^O  +  2k  ^-Kpg(t-t')  ^lr>p  cos(»-e  ) 

J  p2  +  k  (basins' (J>  +  b^cos^ 


(3.19) 


Expansion  of  the  denominator  of  (3.I9)  In  Foiprier  cosine  series  yields 

'3  /;/■'■*  t 

O  m  =  0 

.ir.p  cos  (<t>  -6)  .  .  ... 

i  cos  2m<t>  pdil>dp  , 


(3.20a) 


where 


^0  „  J  p^  +  ik  (b^sin24)  + 


o  ('a  1,  2  _ cos  gm'l’dil’ 

J"  »P  )  “  kJq  p2  +  ik  (b^sinSb  + 


bgcossif) , 


bgCos^f , 


(j.apb) 


Evaluation  of  the  definite  integrals  in  (j.POb)  is  carried  out  in  the 
Appendix  where  it  is  found  that 


»0  <± Tk  - 


\  (± 


,in  .  „-in 


m  -  1 


I 


m-2k-l 


k  =  1 


(3.21a) 


Ski  (±  °  f 


for  m  S  1,  Xi  \2  , 
for  m  S  1,  \i  =  X.2  » 


v;ith  the  notation 
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(5.21b) 


Ni  “  (p®  ±  ikbi)2  , 


^.2  "  (P^  ±  ikbs)®  , 


(kl  -  \£)(\l  +  \2) 


-1 


Recalling  [Wataon  (194k)]  that  the  Bessel  function  of  first  kind  may  be 
represented  by  the  Integral 


T  /  N  1  -Irpcoa*!!  ^ 

Jjj(np)  »  J  S,  ^  cosniiditi  , 


and  taking  41^  ■  0,  (3,20a)  reduces  to 


m  ■  0 


,-/<:p2(  t-t')  T  /  >  j 

i  Jsm  (rp)  Pdp 


(3.22) 


which  appears  suited  to  numerical  evaluation.  At  r=0,  by  (3.20)  and 
with  p'  “  ,  (3.22)  becomes 


(p'+  ikb  )2  (p '+  ikbg)^ 


(3.23) 


which  Is  In  the  form  of  a  Laplace  transform  and  will  be  evaluated  for  several 
specific  shell  geometries  of  the  class  (3.T).  For  an  unlimited  spherical 
shell^  (of  radius  R) 

bi  “  bg  =  -  ,  (3.24) 


and  by  (3.23)  and  a  known  Laplace  transform  [Erdelyi  (1954)]  we  have  at 
r  =  0 


+  Ig  =,  2-n  [-Ci(T)  cos  T  -  Si(T)  sin  ^  |  sin  t]  , 

Ii  -  Ig  -  2«i  [Ci(T)  sin  t.  -  S1(t)  cos  t  +  ^  cos  t]  ,  (3.25) 

T.  -f-  (t-f)  . 

where  S1(t)  and  C1(t)  are  the  sine  and  cosine  Integrals  defined  as  usual  by 


The  unlimited  shallow  spherical  shell  Is  actually  an  unlimited  parabalold  of 
revolution  which  In  the  neighborhood  of  the  apex  may  be  regarded  as  equi¬ 
valent  to  a  spherical  shell  with  edge  sufficiently  far  removed  from  the 
shallow  region  of  Interest. 
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Si(T.)  » y  ^ 


sin  V 


dy  , 


Ci(T:)  »  -  y' 


COB  V 


dy. 


(3.26) 


Also,  for  the  shallow  spherical  shell,  (5. 21)  and  (3.22)  result  in 


/ 


00  -K:p'(t-t') 


+  I2  "  2jr 

_  _  2ffik 

Ii  -  -  — 


■Jn(r.ro')  o'  dp' 


P '®+  ”2 
«  -Kp'(t-t') 


k2 


/o' 


Joitlp')  dp' 


(5.27) 


and  the  solution  is  independent  of  9,  i.e.,  it  is  axiBy..Traetric .  For  an 
unlimited  shallow  circular  cylindrical  shell*  (of  radius  R'  with  generators 
Xg  =  const.} 


bi  -  0  ,  ’’a  -  -  rT 

and  by  a  known  Laplace  transform  [Erdelyi 


!  I 


Kq(+  it')»  t'  =  Ii, 


(I95I+)].  (3.23)  is 

(t-f), 


(3.28) 


(3.29a) 


hence  at  r  «  0 
Ii  +  I2 
Ii  -  I2 

where  K^,  J^, 

( 1944) .  For 


=  Jq(t')  sin  T.'  -  Y^(t.')  cos  T.' 

*  n2  Jq(t')  cos  t'  +  Yq(t')  sin  t' 


(3.29b) 


and  Yq  are  Bessel  functions  in  the  usual  notation  of  Watson 
an  unlimited  shallow  hyperbolic  parabaloidal  shell  with 


=  -  bg  =  b  >  0  , 

by  a  known  Laplace  transform  [Erdllyl  (1954)]  ,  (3.23)  yields 
Ii  =  I2  =  f  [  Hq(t")  -  Yq(t.'')  ]  , 


(3.30) 


(3.31) 


t”  =  kKb(t-t'), 

where  is  Struve's  function  of  first  order  as  defined  by  Watson  (1944). 

The  unlimited  shallow  circular  cylindrical  shell  is  actually  an  unlimited 
parabolic  cylindrical  panel,  which  in  the  neighborhood  of  the  xj,-  axis 
may  be  regarded  as  equivalent  to  a  circular  cylindrical  shell  with  edge 
sufficiently  far  removed  from  the  shallow  region  of  interest. 
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The  solution  for  the  unlimited  flat  plate  (b^  ®  bg  =  O)  cannot  be  obtained  as 
a  limiting  case  of  the  foregoing  solution  for  unlimited  shallow  shells,  since 
for  bj^  =  bg  =  0  ,  by  (3.21) 

Bq  =  ^  ,  am  =  0,  m  S  0  ,  (3*52) 

and  the  integral  in  (3.22)  is  divergent.  This  result  is  due  to  the  fact 
that  the  spectral  representation  (3*2)  for  is  not  valid  for  solutions  of 
the  plate  equations.  For  an  unlimited  plate  Green's  functions  can  be  obtained 
with  the  aid  of  the  extended  definition  of  [see  (4.20)],  but  the  Green's 
function  for  the  bending  solution  will  contain  a  Inr  term  and  an  arbitrary  * 
function  of  time.  In  this  respect  the  Green's  function  for  an  unlimited  flat 
plate  differs  from  those  for  unlimited  shallow  shells  where  the  displacements 
vanish  at  r  =  »  and  no  arbitrary  function  of  time  is  present.  Thus,  a  com¬ 
parison  between  quasi-static  Green's  functions  for  unlimited  plates  and 
shallow  shells  is  not  possible,  and  the  Green's  function  for  the  flat  plate 
will  not  be  derived  here.  For  a  finite  plate,  hov.ever,  no  such  difficulties 
arise. 

o  M 

The  time  variation  of^3  (x,t;  x',t')  at  the  source  point  x  ■■  x' 

vZ  Cv 

(r  =  O)  is  shown  in  Fig.  1  for  the  spherical,  cylindrical,  and  hyperbolic 
parabaloldal  shallov  shells.  These  results  are  obtained  by  numerical  evalu¬ 
ation  of  (3.16),  (3.25)j  (3*29)i  and  (3.30»  for  the  values 

I  =  |,  =  hb  =  ^  ,  V  =  0.3  ,  H  =  0  (3.33) 

where  H  =  0  has  been  chosen  since  H/hk  «  1  in  most  problems  of  practical 
Interest.  For  qualitative  comparison  with  these  quasi-static  Green's  functions, 
since  a  dynamic  Green's  function  for  shells  is  not  available.  Fig.  1  shows  the 
corresponding  time  variation  of  the  dynamic  Green's  functlon^^  (x,t;x't') 
at  r  =  0  for  an  unlimited  flat  plate  as  obtained  by  Nordgren  and  Naghdi 
(1962)  which  may  be  written  as 

[a^]„ .  0  ■  r  ii  (e-f)] , 

1  (5-31+) 

for  ^  ^  “  Ci(iT)  -  Y  -  Si(iT). 

As  seen  in  Fig.  1,  the  time  variation  for  the  three  shells  is  similar  and 
possesses  a  logarithmic  singularity  at  t  =  t'  and  decays  exponentially  as 


-24- 


t  -  t'  -♦  00  .  It  should  be  noted  that  the  quast-statlc  Green's  functions  fall 
to  meet  a  zero  initial  condition,  l.e., 

lim  (x,t;x',t')  5^0,  (5.35) 

t-t'  ->  0+ 

for  general  values'^  of  x  and  x'.  However,  in  the  actual  solution  formula 

(2.55)  ’“>£3  (x,  t;x',t')  8M(x',t')  is  Integrated  with  respect  to  t'  and 

if  ®  (x',t')  is  a  piecewise  continuous  function  of  t'  with  (x',t')  «  0 

M  M 

for  t'  <  0,  then  w  will  be  continuous  in  t  and  vanish  at  t  =  0.  On  the  other 
hand,  if  0  (x',t')  is  an  Impulse  (6-functlon) ,  then  the  singularity  in 

(x, t;x',t')  at  t  =  t'  gives  an  unrealistic  result  for  w,  indicating  that 
the  dynamic  solution  should  be  employed.  In  the  dynamic  solution  for  the  flat 
plate, 3  3  (x, t;x, t')  =  0  at  t  =  t'  as  desired.  Mathematically,  it  is  clear 
that  zero  initial  conditions  cannot  be  imposed  on  the  quasi-static  Green's 
function  since  time  operators  (inertia  terms)  are  not  present  in  the  elasto- 
statlc  equations. 

M 

It  should  be  noted  that4f3  (x,  t;x',t')  given  by  (5-l5)  differ  con¬ 
siderably  from  (x,t;x',t'),  since  by  (2.2),  q  «  tu  when  H/hk  «  1. 

jIj  N  M  ^  ^ 

Thusjgn  will  decay  more  slowly  thanks  as  t  -  t'  -»«  and  for  the  spherical 

and  cylindrical  shallow  shells,  by  (3»l3)»  (3-25)>  (5* 29b),  a  damped 

oscillatory  time  variation  of  significant  amplitude  may  be  expected.  For  the 

hyperbolic  parabaloidal  shell  (3>30),  by  (5'15)  (.3.3l)»^3  (5t,t;x',t')  =  0 

at  X  =  X ' , 

3.2  A  Class  oi^  flectangular  Shallow  Shells.  Let  us  consider  the  class  of  shallow 

shells  (3"T)  now  with  0  S  x^  S  a^  and  boundary  conditions  (2.k)  and  (2.52),  i.e. 


w  =  V2  =  0  ,  Mil  =  H»  ,  Nil  =  N*  ,  on  Xi»  o,ai  , 

w  =  vi  =  0  ,  M22  =  H*  ,  N22  =  N*  ,  on  X2  =  0,82  , 


(3.36) 

(3.37) 


By  (2.19)  and  (3*9),  the  homogeneous  form  of  (5.36)  is  met  if 

similar  situation  occurs  in  the  theory  of  elasticity  where,  e.g.,  the 
thermoelastio  Green's  functions  for  the  entire  space  according  to  quasi- 
static  theory  as  given  by  Parkus  (19' 9)  do  not  vanish  as  t  -  t'  -»0. 
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w 


on  X, 


0,a 


l» 


^ 


0 


6 

H 


a 


O.a, 


a 


(5.38) 


(3.39) 


In  view  of  (5.38)  the  domain  of  the  operators  =  ■^^2  appearing  In  (2.1) 

and  (5.8)  may  be  taken  as  the  set  of  all  Integrable  square  functions  which 
vanish  at  x^  ■  t^ave  piecewise  continuous  and  integrable  square 

second  derivatives  with  respect  to  x  .  Since  by  (3.38)  L  N  vanishes  at 

u  p 

x^  ■  OjSq  »  the  operators  in  (3.8)  of  the  form  be  interpreted  as 

(with  the  domain  Just  specified)  acting  on  L  N  provided  that  N  has  piece- 

P 

wise  continuous  and  integrable  square  fourth  partial  derivatives.  Then, 

since  L,  and  Lp  commute,  we  mcy  again  follow  the  method  established  by  Fried- 
i  ^ 

man  (I956)  and  regard  L  and  >—  as  constants  in  solving  for  the  Green's 
function.  The  result  is  interpreted  using  for  the  spectral  representation 
(3.2)  and  for  and  L2  the  spectral  representations 

5(xi  -xi) 


B(x2-  X^) 


CO 

2 

B  ■ 

y 

sin  - 

mnxi 

sin 

mjTXi 

L 

ai 

ai 

m 

- 

1 

00 

y 

sin 

nnxo 

sin 

njtxo 

^2 

L 

82 

82 

n 

s 

1 

(3.it0) 


respectively,  which  are  given  by  Friedman  (I956,  p.  21*9)*  In  this  manner, 
(3.1)  yields  the  temperature  Green's  function 

Ilk 


^  exp  r-K(m'S+  n'2)(t-t') 

(3.W) 


m  =  1 
n  =  1 


where 


sin  m'xi  sin  m'xj,  sin  n'x2  sin  n'xg  , 


m 


nm 

ai 


nn 


(3.1+2) 


and  a  corresponding  expression  for  Gj^  (x,t;x',t')  follows  upon  replacing 
\  (3.1+1-).  A  similar  technique  applied  to  (3. 11)  results  in 
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(I,  .  I.)  , 

pN  „  2Ehjca  fK^Juit-f)  ^  ^  ^ 

An  A_  ^  ® 


where 


'^)  exp'f‘-/<(in'2+n'^)(t-t')1 


m  =  1  (in'^+  n'^)^  +  lk(b£m'®+  bin'®)^] 

n  ■  1  L 


(3.J^4) 


sin  m'xi  sin  m'xi  sin  n'xg  sin  n'x^ 

By  (3.10),  (3.1^1),  and  the  remaining  components^!  of  Green's 

function are 

00 

jIjN  hKa  t-t ' )  V’  1  /  ,  .v 

=  Ills  ^  2.  ;i^ 

in  =  1 
n  =  1 

^(m'^+n'^)^  (n'^-ym'^)  +  b,k^(in'^-f  n'^)(bpm'^-f  bin'^)  ^ ^ 

^  (m'2+  n'^)'^  +  k2(b2tn'®+  b^^n'^) 

-  exp  •|^-»<(m'^+  n'®)(t-t')  sin  m'xi  sin  n'xs  sin  n'x^  , 

and  a  similar  expression  fori^!  (x,t;x',t')  which  is  given  by  (3*^5)»  with 
M  and  N  as  well  as  subscripts  1  and  2  interchanged.  By  the  same  procedure 
we  obtain  the  Green's  functions 

^  2(i-fv)^a  (li  +  1^)  , 

aiSs  ^ 

„  4(l+v>a  ^-/CT,M(t-t')  y  i(cosm'xi-l) 
aiap  m 

^(m'^+  n'^)(bpm'^+  bin'^(n''^-vm'^)  -  bi(m'^+  n'^)^ 

'  (m'®+  n'^)-^  +  k^(b2in'^  +  b^n'^)^  ^ 

-  l|-  exp  ■r-/<(m'^+  n'^)(t-t')  sin  m'xi  ”'*2  '^'’‘2  > 

a  similar  expression  for^!  (x, t;x',t'),  and 


(3.1^5) 


(3.^a) 


(3.  ^h) 


a^aaki 


m  ^-'^^(t-t')  (j,  _  j.)  ^ 


(3.1+6c) 


where  li  and  are  defined  by  (3. 1+4). 


-27- 


As  an  application  of  the  formulas  of  Section  2  and  the  Green 's  functions 
Just  obtained  for  a  class  of  rectangular  shallow  shells,  we  consider  the  edge 
conditions  (j.j?)  with 


(x,t)  -  f(xi)  , 

=  -f(xi)  . 


X2  =»  0,  t  1  0 


=  -f(xi)  ,  Xg  "  ag  ,  t  S  0  , 

-  0  ,  xi  =  0,ai  , 

«  0  ,  t  <  0  , 


(5.1^Ta) 


and  also  take 


"S  *  "5  ■  On  ■  °  ■ 


(j.Wb) 


where  f(x]^)  is  a  continuous  function  of  integrable  square  on  0  S  x^  S  . 
By  (2.16),  (2.53)1  (3*^7),  we  have 


y»,t)  .  if 

V*'')  ■  X  X  (•‘.'iN'.t')]  0 


3fe  -  0 


which  with  the  Green's  functions  (5*^1)  to  (3*^5)  become 

fl  t’l  =  V  ^  ■  ®*P  +  hn'2  +  qj,)t| 


©N  (x.t) 


32  L  - - 

in«l  m*^+4n'^  +  T]^ 

1  ^ 
n  »  1 


2n*fm  sin  m*xi  sin  2n*X2  $ 

.  ^  M  V  1  -  esP-f-K(m’g+  4n’g+  q„)t] 


ui(x,ti)  -  y  I..- 


n  =  1 


r(ni, '  Ifn '  hn '  ^-ym '  ^ 

1  (m'2  +  l^n'2)4  + 


ym'^)  +  bik^(m'^-fhn'^)(bom'^+  bi4n*^) 
+  k^(b2m'^  +  bi4n'^)® 


1  ?-  -JJJ7  (cos  m'xi  -  1)  fin  sin  2n'x2, 


(3.48) 


(3.49a) 


(3.49b) 


a  similar  expression  for  U2(x, t)  and 
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U„(x  t)  =  -  y  1  -  exp  +  .4n'g_jt,iyt| 

^  ^  82  Z.  m'^  i  4n'2  + 


m  »  1 
n  a  1 


'2  +  +  k2  (b2m'2  +  b 


(m'2  +  4n'2)  (b2m’^  +  bi4n'2)  gn'fm  sin  rn'K^  ain  2n'x2  , 


(5.49c) 


where 


2 

—  /  f(’ti)  si*'  m'xj,  dx^ 
ai  J  o 


(5.50) 


By  (2.57),  (5.9),  (3-43)  and  (3.h7)  the  stress  resultants  are  found  to  be 


r  Nii(x,c) 
1  N22( i) 
(•  Nj^2(^»^) 


Jha  Y*  1  -  exp  -;<(tn'2  + 

^2  L  - ' - — 

m=l  m'^+hn'^  +  Ti 

n  =  1 


>*g  +  4n*^  fm 


+  4n'2  +  Tijj)t1 


(m'2  +  40*2)4  +  k2  (bgin'S  +  bi4n'2)2 

r-'2n*)^  sin  m'xi  sin  2nk2 
m'^an*  sin  m'xi  sin  2n'x2  r 
'  m*  4n'2  cos  m'xi  cos  2n'xp  ’ 


(5-51) 


and  similar  expressions  for  M  „  and  0  are  easily  obtained.  From  the  theory 
of  Fourier  series  it  may  be  shown  in  the  usual  manner  [see,  e.g.,  Friedman 
(1956),  p.271]  that  0„  (3.49a)  converges  nonuniformly  to  the  boundary  values 
(3.47aY  To  this  end  we  write  (3.49a)  as 


=  2^  ^  fm  sin  m'Xj^  sin  2n’x2 


m  “  1 
n  ■  1 


(3.52) 


00 

+  y  A  (m, n)  sin  m'xi  sin  2n'x2 


m  =  1 
n  =  1 


where  A  (m,n)  =  9  (^)  as  n  -»  ».  Thus,  the  second  series  converges  uni¬ 
formly  and  is  zero  at  X2  >>0,82  ,  while  the  first  series  converges  to 

(1  -  ^  )  sir,'  m'xi  fm  »  ( 1  -  f^)  f(xi)  ,  0  <  Xg  <  Sg  (3*53) 


which  approaches  f(xj^)  rnd  -f(xi)  as  Xg  approaches  0  and  82,  respectively. 

In  a  similar  manner  the  remaining  homogeneous  boundary  conditions  (3.36)  and 
(3.47)  way  be  shown  to  be  satisfied  by  (3.49)  and  (3.54)  with  the  series 
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converging  uniformly  in  each  case.  Note,  however,  that  on  Xg  »  0  and 
Xg  =  Hg  ,  Nil  converges  nonuniformly  to  the  values  -Eh  a  and 

Ehaf(xi),  respectively, 

4.  niPROVED  THEORY  OF  SHALLOW  SHELLS. 

GREEN'S  FUNCTIONS  FOR  AN  UNLIMITED  SHALLCW  SPHERICAL  SHELL. 


4.1  Basic  Equations.  Formulation  of  Elastostatlc  Problems.  We  recall  from 
the  note  by  Naghdi  (1956)  that  the  effect  of  transverse  shear  deformation 
may  be  Included  in  the  theory  of  shallow  shells  [in  a  manner  similar  to 
Relssner's  (1945)  improved  theory  of  flat  plates]  by  replacing  (2,l8b)  and 
(2.21)  of  the  classical  theory  with 


and 


"“cp  -  Vp  *  Va  ' 


(1..1) 

(lt.2) 


respectively.  The  components  of  p  represent  the  changes  in  slope  of  the 

cc 

normal  to  the  shell  middle  surface,  and  the  tracer  c  =>  1  is  introduced  in 

(4.2)  to  isolate  the  effect  of  transverse  shear  deformation  in  future  mani¬ 
pulations.  Also,  in  some  equations  it  will  be  possible  to  let  c  =  0  and 
obtain  the  corresponding  equation  for  classical  theory. 

In  the  absence  of  tangential  surface  tractions  (p*^  =  0),  (2.17a) 
leads  in  the  usual  manner  to 

ipP  =  €^4  V.  V  F  ,  (4.3) 

where,  by  the  compatability  equation  associated  with  (2.18a),  together 
with  (2,19a)  and  (4.3),  the  Airy  stress  function  F  must  satisfy 

F  +  Ehe^  6^^  V  V_  Z  V  V  w  +  Eh  a  =0  (  4. 4) 

ct  p  X  q  w 

just  as  in  classical  theory.  Introduction  of  (2.19b),  (4.1),  (4.2),  and 

(4.3)  into  (2.17b,c)  results  in 

+  v%)  +  ,  (1..5) 

+  p  =  0, 
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(4.6) 


ie( i+v)c 


O^+V^w)  -D 


2 

+  D  a  V  “  0- 

rl 


which  contain  the  effect  of  transverse  shear  deformation.  Next,  by  the  two- 
dimensional  version  of  the  Stokes-Helmholtz  theorem  of  vector  analysis  we 
may  write 

V  ■  (“-T) 


and  then  (4.6)  becomes 


2  _i 


V  H'  -  X  (H'  +  w)  -  (1  +  v)  a9^  -  f, 


"k.  10  k'  -  a 

^  h^c  ^  1-v  *  ’ 


where 


.  ch® 

•  50^  • 

Letting 

H  -  H'  +  \£  .  K  =  K'  +  Xg, 

by  (4.9)>  equations  (4.7)  and  (4.8)  read 

B  -  V  H  +  a  e’l^  V,  K, 
a  OT)  X  ’ 

V^H  -  x‘^(H+w)  -  (1  +  v)a0„  =  0, 

w 

V  K  -  K  =  0 
h‘='c 


(4,8) 


(4.9) 


(4.10) 


(4.11) 

(4.12) 
(4.15) 


which  are  the  same  for  all  shallow  shells,  and  (4.5)  is 

XD(V^H  +  V^w)  +  ^  V^Vp  Z  F  +  p  -  0.  (4.14) 

Thus,  all  solutions  in  the  improved  theory  of  shallow  shells  may  be  repre¬ 
sented  in  terms  of  F,  w,  H,  and  K  which  are  governed  by  (4.4),  (4.12),  (4.13), 
and  (4.14).  It  may  be  noted  that  for  flat  plates  (Z  =  O),  (4.12)  and  (4.13) 
remain  unchanged  while  (4.4)  reduces  to  an  equation  for  F  alone  and  (4.14) 
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reduces  to  an  equation  containing  H  and  w.  For  bending  of  flat  plates  the 
system  (4.12),  (4.15),  and  (4.14)  with  2  =  0  is  equivalent  to  the  formulation 
by  Reissner  (1945)  which  Involves  a  stress  function  representation  for  Q*^. 


4.2  Green's  Functions  for  an  Unlimited  Shallow  Spherical  Shell  Under  Normal 

Traction  and  Temperature  Field.  For  a  shallow  spherical  shell  (of  radius  R) 
referred  to  rectangular  Cartesian  coordinates  x 


-x^x 

a  a 


(i^.l5) 


the  basic  equations  (4.4)  and  (4.14)  simplify  somewhat,  while  (4.11)  and 
(4.12)  remain  unchanged.  For  an  unlimited  shallow  spherical  shell  the 
solution  may  be  written  in  the  integral  form 

poo  poQ  r 

w(x)  ”  /  J  w  (x;x')  p(x')  +  w^(x;x')  6^  (x') 

*^00  *09 

(4. 

+  w^(x;x')  ®jj(x')  dxi  dx^  , 

and  similar  expressions  for  F,  H,  and  K,  where  the  three  sets  of  Green's 
functions  wp,  Fp,  Hp,  Kp;  and  w^^,  F^^,  Hj^,  are  solutions 

of  the  basic  equations  with  p  =  &(x-x'),  =  &(x-x'), 

p  =  =  0;  and  =  fi(x-x'),  P  =  respectively.  Thus,  all  three 

sets  of  Green's  functions  may  be  obtained  from  solution  of 


(4.16) 


2  2  yu  2 

VV  ■  R  “g  "  ■  N  > 

\(v\  +  v\)  -  ^  ^  i  6(x.x'). 

V  Hq  -  (Hq  +  Wg)  ■  (1  +  v)  6(x-x')  , 

2  10 

V  K-  =  0  , 

G  n  c  G 

where  p,  ,  and  are  constants  taken  in  turn  as  p  «  1,  =  0; 

®  1,  p  =  Sjj  =  0;  and  0^^  =  1,  p  =  0^  ■»  0,  for  each  set  of  Green's 
functions. 

To  solve  (4.1?)  we  Introduce  polar  coordinates  r  and  0  through 


(4.17a) 

(4.17b) 

(4.17c) 

(4.17d) 


Xi  -  s»  r  cos 


=  r  sin  0, 


(4.18) 


-32- 


and  recall  that 


S(x-x')  =  —  8(r)  , 


(iM9) 


and  the  definition  of  the  Laplacian  may  be  extended  to  functions  U(r)  not 
regular  at  r  =  0  according  to 


'5'  ■■  r^O  '  r^O 


(4.20)' 


lim  r  KgJ  =  0 
r  ->0 


(421a) 


where  So/5r  indicates  differentiation  in  the  ordinary  (not  symbolic)  sense. 
With  the  aid  of  (4.19)  and  (4.20),  (4.17)  may  be  written  as 

[/•?%  -  f  ^  [u«  r  (T%  -  f  »^)] 

+  I" ^  f  ^ 

“-r  0  -*  *-r  ->  0  ° 

+  ^  [u.  rF^.O  , 

i  v\]  +  ^  [  11. Jj  (>.-7h,  +  w^) 

(4.2 

■  RD  r  ^  '^G^  '  RD  ^G^  " 

'-r  -*  0  -> 

[y\  -  x.'‘  (h^  +  "j)]  +  ^  [u.  t  ^ 

h'  “a]  -»■ 

•J  '-r  -» 0 

C'^’^'k  -  K  1  +  r lim  r 

I  G  h^c  Gj  r  Lr_o 

r  -I 

+  SJil.1  lim  r  K.  =  0 
r  L  Gj 


(4.21b) 


(4.21c) 


(4.21d) 


which  are  satisfied  if  each  expression  in  brackets  vanishes.  The  first 
brackets  in  ( 4.2la,b, c,d)  lead  to  the  general  symmetric  solution 

*-X-  "" 

The  corresponding  extended  definition  of  the  three-dimensional  Laplacian 
Is  derived  by  Friedman  (1959)>  P"  255. 


(4,22a) 


Wg  =  ~  Bi  -  Bg  Inkr  +  “  1)  B3I0  (  ^ir) 

+  B4K0  (^ir)  +  -  1)  Bsl„  (5gr)  +  Bg.lC„(53r) 

Hg  =  Bi  +  Bs  Inkr  +  B3  I,  ( +  B4K„(?ir) 

+  B5l„  (?E  r)  +  BgK„  (^sr)  , 

Pg  -  ^  +  Inkr  +  RD  Us  5i  1.  (£ir)  '  ^  ' 

+  84^1  Ko(?ir)  +  lo(^Er)  +  BgK„  ( . 

where  Bi,  Bg,  .  .  .  B  are  constants  and 


(4.22b) 


( 4.22c) 


(4.22d) 


R.£.,  £=>.0, 


(4.25) 


Recalling  that 


lo(Cr) 


,  as  r  ->  CO 


Re  ^  >  0  , 


(4.24) 


then  In  order  that  w.  and  K_  have  at  most  a  singularity  of  0( Inr)  as 

u  Ca 

r  -»  00  ,  B3  =  B5  =  B^  =  0,  and  without  loss  in  generality,  B^  *»  B^  >=  0. 
Further,  since 

lim  r  Ko(5r)  =  lim  Ko(5r')  =  0, 

r  -»  0  r  -)  0 

lim  r  Ko(5r)  =  -1  ,  lim  r  /  K,(£r)  =  -  , 


(4,25) 


r  0 


(4.22)  meet  (4.21),  provided 


r  -»  0 


■  2.(£;4  £i)  [d  *  ’“'““‘’n  +  E=  <1  +  v)  a?J  . 

“6  ■  Mil  -  £il  [d  ■"  ““H  ■* 


(4.26) 


Thus,  the  solution  (4.22)  which  determines  the  Green's  functions  is 


G  m 


J  +  Rk  aGjj  +  ^2  ( 1  +v)  aQ^ 

|_  §  +  Rka\  +  (1  +  v)  cc\\ 

pR  +  Eh  (4.27) 

1  2  - 

+  RD^s  (1+v)  a\J  Ko(tir)  -  ^1)  [pR  +  Eh  +  RD(;i(  1+v)q@j^ 


,  {Kl-  .U 


Ko(?ir) 

K»(Ur)  5 


""g  “  2«r5T”TiT 


(£  +  RKae^  +  (I4v)  a0^) 


-  (  J  +  Rk'^aSjj  +  (1+v)  a^j^)  K<,(^2r) 


To  obtain  corresponding  results  for  the  classical  theory  of  shallow  shells  we 
neglect  the  effect  of  transverse  shear  deformation  by  setting  c  =  0,  in  which 
case  for  ~  ^  (^.P?)  reduces  to  the  solution  given  by  Reissner  (1946) 

for  an  unlimited  shallow  spherical  shell  under  concentrated  load,  and  for 
p  =  0  (4.27)  reduces  to  the  solutions  given  by  FlUgge  and  Conrad  ( I956)  for 
concentrated  temperature  resultants. 

With  the  aid  of  (4.25),  in  the  neighborhood  of  r  =  0  (4.27)  is 

Wq  =  ^  (  ^  +  Rk^a  )  |^(l-X.  k  )  r^  In  kr  -  4  In  kr  +  . . 


.  (1+^)  Q  B 

2Jt  M 


1242 

In  kr  -  -jj-  X  k  r  In  kr  +  . . . 


"g“  - 


4 

R\k  - 

—  P 


■  2 

r  In  kr  +  . . . 

Eh  a0jj 

•« 

2”  L 

1  22 

+  -jjl  \k  r  In  kr  + 


2 

.  (1+v)  aRDk  -S 

•••  - 


In  kr 


r  In  kr  + 


(4.28) 


“0  ■  •  8^  '  D  + 


In  kr  +  . 


(1+v)  g  e 


2jt 


M 


In  kr  +  . . , 


which  differs  from  the  classical  theory  (c  =  X  =  O).  In  particular,  for 
the  case  ®  according  to  the  Improved  theory  w^  -»  «  as  r  -»  0, 

while  in  the  classical  theory  Wp  -»  0  as  r  -»  0.  In  this  connection  we 
recall  that  in  most  other  derivations  of  singular  solutions  for  concentrated 
loads  in  the  classical  theory  of  plates  and  shells,  the  requirement  Wp  ->  0 
as  r  0  is  imposed  a  priori;  however,  no  such  requirement  is  needed  in  the 


Ko(Ur), 
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foregoing  treatment.  Indeed,  the  requirement  Wp  0  as  r  ->  0  for  the 
Improved  theory  would  lead  to  an  incorrect  result  for  the  Green’s  function 
as  seen  from  (4.28).  It  may  be  recalled  that  Green's  functions  also  are 
unbounded  at  the  source  point  in  the  three-dimensional  theory  of  elasticity. 

To  see  that  the  p  term  in  (4.27)  represents  a  concentrated  normal  force 
at  r  =  0,  by  equilibrium  considerations  the  total  resultant  normal  force 

Q 

V“.  on  a  vanishingly  small  circular  region  with  center  at  r  =  0  is 

=  -  111..  r  (Q®  +  fN^)rd0.  (4.29) 

r-^0"^0 


By  (4.2),  (4.5),  and  (4.11)  in  polar  coordinates  with  polar  syntnetry,  we  have 

(4.30) 


s<>'g+”g)  • 


N®  = 

r  Sr 


and  by  (4.28),  (4.29)  yields 


V 


G 

N 


=  P 


(4.31) 


as  desired. 

Green's  functions  for  a  flat  plate  with  the  effect  of  transverse  shear 
deformation  Included  cannot  be  obtained  as  a  limiting  case  of  (4,27)  as 
R  -»  «  .  Instead,  we  must  return  to  (4.21),  let  ^  =>  0  ,  and  then  in  the 
same  manner  as  for  the  shallow  spherical  shell  we  obtain 


In  ^  -  4\  (1  +  In  ^  ) 


( 1  +  v)  g 

2n 


0 


M 


«G 


r 


+ 


( 1  +  v)  g 

2ji 


M 


(4.32) 


Eh  g  ,  r 
— —  y  In  — 

2jt  N  h 


When  the  effect  of  transverse  shear  deformation  is  neglected  (c  =  X  =  O) , 

(4.32)  reduces  to  known  results  for  classical  plate  theory.  For  the  flat 
plate  the  Improved  and  classical  theories  give  different  fundamental 
singularities  only  for  the  Green's  function  w  . 
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5-  APPENDIX: 

EVALUATION  OF  CERTAIN  DEFINITE  INTEGRALS 

Here  we  shall  evaluate  the  definite  Integrals  (5.20b)  by  contour  integration 
around  the  unit  circle  in  the  complex  Z  plane  where 

Z  =:  r  i  .  (5.1 

Thus,  on  the  unit  circle  C  (r  =  1) 


sin  $  = 


JL  i. 

Z2  ,  Z  2 


cos  '!>  = 


1  -i 

H2  +  2  2 


cos  2m<!>  « 


Z™  +  z'"* 


dZ  =  2iZd(li  , 


and  for  b^^  4  t>2  (3*20b)  becomes 

“  ^i(xi  -  Xf) 

(±  A 


"m  (± 


(Z  -  p)  (Z  - 

(Z*"  +  z'"^)  dZ 
(z  -  p)  (z  -  3-1) 


2ii. 

ez. 

[  + 

\p 

1  -  XiXp  ■  XiX 

>^1 

+ 

Xg 

+  ^j^X 

where  \i,  \2»  3  defined  in  (5.21b).  Consider 


|P|^  =  pp  = 


where  p  denotes  the  complex  conjugate  of  p.  We  may  write  and  \g  in  the 
polar  form 

,  2,  A.i'i  ,  2.  _  ,jt 

\i(+  k,p  )  =  Rj^i”  ^  ,  X2  (+  k,p  )  “  R2  ^  , 


-  ^  ®  ^  i  > 


whence 


^1^2  ^  (^1  "  ^2)  I 

-  2  s  ($1  -  V?)  s  I  . 

Therefore,  by  (5'*+)  bo  (5»6) 

\iXp  +  ^  0  , 


'el  s  1  , 


This  method  is  used;  e.g.,  by  Churchill  (19^) 
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where  the  equality  signs  apply  only  when  p  =  0  for  a  special  subclass  of  (5*7) 
Excluding  for  now  the  case  p  =  0  ,  the  integrand  in  (5.2)  has  a  pole 
inside  the  unit  circle  only  at  g  =  p,  hence  by  the  residue  theoran  we  have 


Similarly,  with  the  aid  of  the  expansion 


-  p)  (g  -  p-i)  ~  z  z  ^ 


J-2k  J 


2"  , 


j  =  0  k  =  0 


(5.3)  becomes 


2  A™  +  a""*  li 

•«  -  Ki) 


m  -  1 

I 

k  =  0 


P 


m-2k-l 


m  g  1  ,  bj^  5^  bj 


(5.8) 


(5.9) 


(5.10) 


For  the  special  case  p  =  0,  since  (+  k,  p^)  appears  in  the  integrand  of 
(3.2),  it  is  sufficient  to  take 

%  k»P^)  > 

p  0 

where,  by  (5.9)  and  (5.I0)  the  limit  in  (5. 11)  exists  for  b^^  /  bg.  For  the 
spherical  shell  b^  =  bg,  by  (3.21b),  =  \2  and  (3.20b)  yields  (5-8)  and 

(±  k,P^)  =0,  m  S  1.  (9'12) 

Thus  we  have  obtained  the  results  expressed  by  (3.2la). 
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Figure  1.  Time  variation  of  the  Green's  £unctlonQr3  (xit;x';t')  at 

X.  ■  x'  (r  =  O)  according  to  quasi-static  theory  for  unlim^tedj^spherlcal^ 

cylindrical,  and  hyperbolic  parabaloidal  shallow  shells  (  ^  ^  « 

V  =  0.5)  and  for  qualitative  comparison  according  to  dynamic  theory 

1  )  for  an  infinite  flat  plate. 
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